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Abstract 

I present a Ginzburg-Landau theory for hexagonal oscillations of the upper critical field 
of UPt3 near T c . The model is based on a 2D representation for the superconducting order 
parameter, fj = (771,772), coupled to an in-plane AFM order parameter, rh 3 . Hexagonal 
anisotropy of H c 2 arises from the weak in-plane anisotropy energy of the AFM state and 
the coupling of the superconducting order parameter to the staggered field. The model 
explains the important features of the observed hexagonal anisotropy [N. Keller, et al, 
Phys. Rev. Lett. 73, 2364 (1994).] including: (i) the small magnitude, (ii) persistence of 
the oscillations for T — > T c , and (iii) the change in sign of the oscillations for T > T* and 
T < T* (the temperature at the tetracritical point). I also show that there is a low-field 
crossover (observable only very near T c ) below which the oscillations should vanish. 



I. INTRODUCTION 

The heavy Fermion superconductor UPt3 has attracted a great deal of attention because of its remarkable 
low temperature phase diagram. The key features are: (i) the existence of two superconducting phases in zero 
field with a small separation of the transition temperatures, AT C ~ 50 mK , compared to T c ~ 500 mK , (ii) 
the existence of three superconducting phases in a magnetic field (iii) with phase boundaries that meet at a 
tetracritical point (T*, if*) on the upper critical field line. |l],^] Interpretations of this phase diagram fall into 
two main categories of theoretical models: (i) theories based on two symmetry un-related order parameters 
which are accidentally nearly degenerate, gj4|] and (ii) models based on a single multi-component order 
parameter belonging to a higher dimensional representation in which the degeneracy is lifted by a weak 
symmetry breaking field (SBF). |3]-[l0| UPt3 is a close-packed hexagonal crystal with two formula units per 
unit cell; the point group is D^h- Two candidates for an intrinsic SBF have been identified experimentally. 
Neutron scattering experiments show that UPt3 develops an anti-ferromagnetic order parameter, m s _L c, 
below a Neel temperature of Try — 5 K. [jllj The magnitude of the ordered moment is small, ~ 0.02/jb per U 
atom, and is directed in the basal plane thus breaking the in-plane hexagonal symmetry. Evidence in support 
of an AFM symmetry breaking field was found from pressure studies of the superconductivity and AFM order. 
Heat capacity measurements show that the splitting of the superconducting transitions is suppressed under 
a hydrostatic pressure of p c ~ 3.8 kbar, fl2| while neutron scattering measurements show the supression of 
AFM order at roughly the same pressure, p* ~ 3.2 kbar. Jl3| However, the observed magnetic Bragg peaks 
indicate finite-range AFM order with a correlation length £ a /m — 250 A depending on the crystal, which 
has led to discussion of whether or not the AFM order is intrinsic property of UPt3. jl4j Alternatively, 
a macroscopic strain field was proposed jL5| as a possible SBF based on the experimental observation of 
a complex incommensurate structural modulation in UPt3 with characteristic wavelengths of order 10-50 
lattice spacings and typical domain sizes of order 10 4 A. However, these structural modulations seem unlikely 
to be the principal SBF responsible for the double transition since the high formation temperature suggests 
that they persist to pressures well above that required to destroy the double transition and the AFM order. 

In a recent article, Keller, et al. |17| reported the observation of a weak hexagonal modulation of the upper 
critical field of UPt3 as a function of the orientation of the field in the basal plane. The key features of these 
measurements are: (i) the small magnitude of the oscillations, \5H^2 X / H C 2\ < 0.02, (ii) persistence of the 
oscillations for T — > T c , and (iii) the change in sign of the oscillations for T > T* and T < T*, where T* 
is the temperature at the tetracritical point. In this paper I show that hexagonal oscillations such as these 
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are a consequence of the in-plane AFM symmetry breaking field above T c that is responsible for the double 
transition in UPt3. 

II. PERFECT HEXAGONAL SYMMETRY 

In order to sharpen the argument I briefly review the important theoretical results on the anisotropy 
of the in-plane upper critical field in any perfectly hexagonal superconductor. First consider a hexagonal 
superconductor described by a one-component order parameter, 7/. This case includes conventional super- 
conductors for which r\ is invariant under D$h, as well as unconventional cases where 77 changes sign under 
one or more symmetry operations. In the Ginzburg-Landau (GL) limit the free energy functional for any 
superconducting order parameter belonging to a one-dimensional representation is 

T = a{T)\n\ 2 + /3|r;| 4 + K||(A»?)(A»?)* + K ± (D zV )(D zV y , (1) 

where D% = V, + 2e/hcAi and repeated indices are summed over (x,y). An important point here is that 
the GL functional has an accidental symmetry; to second order in the gradients there are no terms in the 
functional that differentiate hexagonal symmetry (D^h) from cylindrical symmetry (Z?oo/i)- Consequently, 
the upper critical field is independent of the orientation of H in the basal plane; = |||a(T)|/^//«[j7ti ~ 
(1 — T/T c ). Hexagonal anisotropy of H^r 2 shows up only at lower temperature when higher-order gradients 
become significant. The lowest-order contribution to the free energy functional that is invariant under D§h, 
but not Dooh, is 

Fhex = K 6 \(D x + iD v fr]\ 2 . (2) 

It is straight-forward to show that this term leads to hexagonal oscillations of H^ 2 (T) of the form, 

8H^(&, T) oc « 6 (1 - T/T c f * cos(6t5) , (3) 

to leading order in Kq. The key point is that the hexagonal anisotropy vanishes as (I — T/T c ) 3 as T — > T c . 
This result is not limited to superconductors with a one-dimensional order parameter. 

Burlachkov |Q showed that H^r 2 is isotropic for fields in the basal plane, in the GL limit, for any of the 
two-dimensional representations of Dgh- In this case the GL functional can be written as 

Tsc = a{T)\ ni \ 2 + Al^l 4 + f3 2 \n 2 \ 2 + k a {D zVl ){D zJli )* 

+ Kl (D i r 1j )(D i r, j )* + K 2 (D i r h )(D j ri j )* + k 3 (A^)(A%)* > [ ' 

where the order parameter (771, rj 2 ) transforms according to one of the 2D irreducible representations (Ei or 
E 2 ) of D 6 h. The isotropy of H^ 2 in the GL limit again follows because to second-order in the gradients and in 
(r]i,r]2) the GL functional is invariant under the larger group, -Dooh- Again, one has to examine sixth-order 
terms (therefore, of order (1 — T/T c ) 3 ) in order to develop hexagonal anisotropy in H^. JTofl 

Recently, Mineev examined GL models based on two accidentally nearly degenerate ID representations. 
Such a models have been investigated by several authors as candidate theories of the H-T phase 
diagram of UPt3. An important feature of these models is that they do not rely on the coupling of the 
order parameter to a SBF that reflects a weak breaking of hexagonal symmetry above T c . The basic result 
for the anisotropy of is the same; hexagonal oscillations vanish as (I — T/T c ) 3 near T c . |2(J Thus, the 
experimental observation of hexagonal oscillations of in UPt3 for T — > T c is in conflict with any GL 
model of the superconducting phases that is based on perfect hexagonal symmetry above T c . 

Below I show that hexagonal oscillations of H^(T) do appear in the GL limit (T w T c ) for the class of 
models based on a 2D superconducting order parameter coupled to an AFM symmetry breaking field present 
above T c . It may initially seem odd that one obtains six- fold oscillations from an AFM order parameter 
that reduces the symmetry from hexagonal to orthorhombic. Six-fold oscillations of are shown to be a 
natural consequence of a small in-plane anisotropy energy of the AFM order parameter, while orthorhombic 
anisotropy of H c2 appears in the limit of large in-plane anisotropy. 
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III. BROKEN SYMMETRY MODEL 



The free energy functional for a 2D superconducting order parameter, ff, coupled to an AFM order param- 
eter, m s , is the sum of (i) T sc , which represents the GL functional for the superconducting order parameter, 
(ii) Tajm-, which is the free energy for the AFM phase above T c , and (iii) !F sc - a f m , which represents the 
coupling between superconductivity and anti-ferromagnetism. The general form for the superconducting GL 
functional in the absence of a SBF has been discussed by many authors pi]] , and is given in eq. [|. First, 
consider the magnetic free energy above T c . 



A. Free Energy of the AFM state above T c 



For T < Tn the mean-field approximation for the magnetic free energy functional should be reliable 
provided the correct magnetic ordering has been identified. I assume an order parameter, m s = (m x , tn v , m z ), 
describing AFM correlations between neighboring U atoms in a double unit cell, and start from a Landau 
functional for the magnetic free energy density, 

Faf m = a(T)\m s \ 2 +b\rh s \ i + ... _ 

+a z m 2 z + d Re(m x + im y ) 6 + g(m s ■ H) 2 . 

The first line of terms is invariant under the full spin-rotation group. They are primarily due to exchange in- 
teractions in typical magnetic materials. Thus, neglecting the anisotropy terms, one has mo — y/\a(T)\/2b oc 
\T-T N \^ 2 for the magnitude of the AFM order parameter in the exchange approximation. 

The third term represents the leading order uniaxial anisotropy energy. Anisotropy energies arise from 
spin-orbit interactions and are typically small compared to the exchange terms. I assume a z > which 
favors in-plane AFM order. 

The sixth-order term in m s is the leading term in a GL expansion for the in-plane anistropy energy. The 
in-plane anisotropy energy is assumed to be a small perturbation to the AFM exchange energy, i.e. 

Uanis =dm% < U exch = -\a(T)\ml . (6) 

In the absence of an external field the in-plane anistropy energy leads to six degenerate minimum energy 
orientations for rh s , which are equivalent to three commensurate wavevectors in a spin-density wave descrip- 
tion. There are two sets of preferred orientations depending on the sign of the coefficient c'. If c' > the 
moments prefer the set of alignments, {i?„ = nn/3 + ir/6;n = 0, ...,5} that includes the a* axis, while for 
d < the set of minimum energy orientations are shifted by 7r/6, and includes the a axis [I use the notation 
inRxjf. @]. 

The Zeeman energy is quadratic in H for an antiferromagnet and prefers the AFM order parameter to 
be aligned perpendicular to the field (g > 0). Thus, for a general orientation of the field in the basal 
plane there is competition between the anisotropy energy and the Zeeman energy. The field at which the 
anisotropy energy is comparable to the Zeeman energy defines a cross-over field scale, 




At low fields, H <C H an i S , the AFM order parameter is essentially locked by the in-plane anisotropy energy to 
an equilibrium orientation, i? n . However, for H > H an i S the Zeeman energy dominates the anisotropy energy 
and the orientation of the magnetization will adjust to remain approximately perpendicular to H . This case 
is particularly relevant to the discussion of the hexagonal oscillations of H C 2- In the limit H H an i S the 
orientation of fh 8 is fixed by the field, •& = dn — tt/2, or m s — m s (sini?ij, — cosflu). The magnitude of the 
staggered moment is then to a very good approximation determined by minimizing the AFM free energy 
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T afm = a(T) to 2 +bm i s - c' m 6 s cos(6i? H ) , (8) 

at fixed orientation $ = du ~~ 7r/2. Treating the anisotropy energy as a perturbation leads to a small correction 
to the exchange approximation for the AFM order parameter; m s = too + Sm s . Retaining the leading order 
corrections to the stationarity condition gives a hexagonal modulation of the AFM order parameter, 



Sm s 3 / U a 



too 4 V U, 



J exch 



cos(6$ H ) , H > H ams . (9) 



This modulation of the AFM order parameter by the anisotropy potential leads to a hexagonal modulation 
of the upper critical field in the SBF model for the double superconducting transition. 



B. Coupling between Superconductivity and AFM 

I assume a specific 2D model, based on an E2 order parameter, which I discussed recently in the context 
of the H-T phase diagram, pi In this model the gradient coefficients n 2 = K3 ~ for weak hexagonal 
anisotropy]^] however, for in-plane magnetic fields the principal arguments that follow apply to any of the 
2D models. The parameters of the GL functional in eq. ^ are calculated from Fermi-liquid theory. [ p3|]24| ] 
The gradient coefficients, k,\ ~ iV/(w^) 2 /(7rT c ) 2 and K4 ~ Nf(v^) 2 /(TtT c ) 2 , are determined by the density 
of states, Nf, and the anisotropic components of the Fermi velocity, '". These coefficients determine the 
uniaxial anisotropy of H c \ and H C 2- The other coefficients, a(T) = Nf(T — T c ) and 02 = \Pi — 16^2 -pf , 
determine the condensation energy; T c is the transition temperature in absence of coupling to the AFM 
order parameter. 

To calculate the hexagonal anisotropy of H C 2 near T c I need the the terms in the free energy that describe 
the coupling of the superconductivity to the AFM order parameter. The leading order symmetry breaking 
terms are 

Fsc-afrn - -6 TO 2 NfT c (\ m | 2 - \ m \ 2 ) 



2 «i(|X>x»7i| a -|-Dx»72| a ) ( 10 ) 



-ej_ to,- 

+e N m 2 n 4 (\D zm \ 2 ~\D zm \ 2 ). 

The first term is responsible for the double transition in zero field, while the gradient terms generate an 
asymmetry in the slopes of the two branches of H C 2 for T > T* and T <T*. 

Combining the coupling terms with the superconducting free energy (eq. |J) gives, 

J . sc ~t~ J. sc—afm ^+ 

+4(D ±Vl ) ■ (D±r]i)* + Ki0^ 2 ) ■ (D ±m )* (11) 
+4{D z r h )(D z r h )* + ^{D zm ){D zm )* , 

with a±(T) = Nf(T — T^). The fourth-order terms are omitted since they are not relevant to the determi- 
nation of H C 2- The SBF field generates a splitting of T c and the gradient coefficients, 

T± = T c (1 ± £to 2 ) ; Kf 4 = K 1A {l±e ul m 2 s ). (12) 

The upper critical field is obtained from solutions of the linearized GL equations, 



1 In this context the term 'weak hexagonal anisotropy' refers to the Fermi surface. This anisotropy is the origin of 
the hexagonal anisotropy of H^ 2 {T) in perfect hexagonal systems, which vanishes for T — > T c as (1 — T /Tc) 3 ■ This 
same Fermi surface anisotropy is responsible for the deviation of K2,3 from zero in the E2 representations. fllCl] 
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K l D ± V2 + K4 D\ T] 2 = a- (T) 772 . 
The resulting value of H^(T) is then the larger value of the solutions for the two branches, 

hcN f (T±(d H )-T) 



2c 



(13) 



(14) 



(tetracritical point). The 
denotes the branch with T < T*. The two branches 



For Ki k~1 > k^k^ the upper critical field has a kink at the temperature T 
solution H~c2 denotes the branch with T > T* and H~ 2 
correspond to transitions from the normal state to a superconducting state with order parameters, ff ~ (1, 0) 
for T > T* and ff ~ (0, 1) for T < T*. Note that Hf 2 (T) depends implicitly on the orientation of the field, 
§Hi through the AFM order parameter m s . 

Hexagonal oscillations of (T) result from the modulation of the AFM order parameter by the sixth-order 
anisotropy energy. In the limit H 3> H an i Sl the extrapolation points, T c ± (^^f ), show a hexagonal modulation. 
The modulation can be scaled in units of the zero-field splitting of the transition, AT C /T C — e (\ + (3\ / fa) : 
and the ratio of the anisotropy energy to the exchange energy, 



T±{$ H )=T C 



AT r 



U„ 



±3 

2 l + Pi/fc \U exch 



[cos(6tf ff ) - 1] 



(15) 



The parameters (without arguments) denote the exptrapolation points in zero field. Note that the 
modulation of the two extrapolation points, T+ and T~, is out of phase by 180°. This explains the change 
in phase of the oscillations of H^(T) for temperatures above and below the tetracritical point. Figure 1 
shows the upper critical field for two orientations of the magnetic field. Note the change in H^(T) upon 
rotating the field from = (minimum of the anisotropy energy) to $h — tt/6 (maximum of the anisotropy 
energy). This change in phase of the oscillations for T > T* and T < T* has the same origin as the splitting 
of T c in zero field by the AFM symmetry breaking field. Q 
The oscillations are small in magnitude, 



1^ 



= i \H±{fi H = 0) - H%(6 H = 7T/6)! 



dT 



3 

* — 



2 \Uexch 



U„ 



(16) 



and weakly dependent on temperature, except for T « T* , where the oscillations change sign (see Fig. 
1), and for temperatures very close to T c where H C 2(T) drops below the anisotropy field, H a ms- From the 
magnitude of the hexagonal anisotropy measured by Keller, et al. [ p5[ , \6H C 2\ max ~ 10 -3 T for T ~ 507 mK > 
T* ~ 430 mK, one obtains an estimate for the ratio of the anisotropy energy to the exchange energy. Taking 
AT C ~ 60 mK, dH+/dT ~ 4.5 x 10~ 3 T/mK, @ gives U anis /U exch ~ 0.015. 

In addition to the modulation of the extrapolation temperatures, the slopes of the two branches of H C 2(T), 
are modulated by the anisotropy potential, 



S ± (^ H ) = -dHfJdT = 5 1 T (e,| + e ± ) 



3 fU a 



4 V U, 



exch 



[cos(6i?jj) - 1; 



(17) 



where S is the slope corresponding to a minimum of the anisotropy energy. These terms generate temper- 
ature dependent amplitudes for the oscillations of the two branches of H^(T) as shown in Fig. 1. 



2 Note that Fig. 1 includes the additional corrections discussed in the next section; however, the dominant effect is 
the modulation of the extrapolation temperatures by the anisotropy energy. 
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C. Higher-order corrections at lower fields 



The oscillations of H^(T) obtained above are for fields large compared to the anisotropy field, H anis . 
While it is reasonable to assume that H an i S is very low, it is non-zero and so provides a cross-over field 
below which the hexagonal oscillations of should vanish. For H > H an i S there are also corrections to the 
'high-field' limit discussed above. These corrections, as well as the low-field cross-over behavior, are obtained 
from the minimization of the AFM free energy functional with respect to both m s and d for fixed Sh- The 
stationarity condition 5T a f m / 5m s = generates, 

I f H 1 




cos(6i?) + - cos^tfn - 0) } , (18) 



m 4 \U exch J 3 V H, 



to leading order in U an i S /U exc h. The second term in eq. [l^ originates from the Zeeman energy. The relation 
between the orientation of the applied field (du) and that of the AFM order parameter ($) is determined 
by the stationarity condition 5!F a f m / 5§ = 0, 

sin(6tf) = \ (j^—) sm(2(ti H - 0)) . (19) 

In the zero-field limit one obtains the expected result that the AFM order parameter is locked to the lattice 
at the minima of the anisotropy energy given by nir/3 + ir/6 for c' > 0; the solutions wr/3 correspond to 
maxima of the anisotropy energy. In the opposite limit, H an i S /H — > 0, the AFM order parameter is oriented 
at ±7r/2 relative to the field; the other solutions at ±7r corrrespond to maxima of the Zeeman energy. 

For H anis /H <C 1 the orientation of m s is no longer strictly orthogonal to H. The deviation, (3 = 
(dn — "&) — tt/2, oscillates with an amplitude proportional to (H an i S /H) 2 , 

/3^3(^p) sm(6V>. (20) 

This phase shift vanishes near a minimum of the anisotropy energy and is maximum at orientations of H 
corresponding to a maximum in the slope of the anisotropy energy (i.e. maximum anisotropy torque on rh s ). 
This correction, combined with the correction from the Zeeman energy, gives 



Sm s _ 3 / U a 
too 4 \U e . 



(21) 



TT \2 \ / 77 \2 



f(0 B ) = cos \6#b 18 sin(60ff) I - 3 j si^(6^ H ) , (22) 

for the modulation of the AFM order parameter parameter. The corrections from the Zeeman energy and 
the angular deviation from tt/2 generate non-sinusoidal hexagonal oscillations of H C 2(T). Figure 2 shows 
these oscillations for temperatures above and below T* . Note the sharpening of the oscillations for H\\a* 
({>rj = 7r/6). A similar effect is also seen in the experimental results for H^T) reported by Keller, et al. 

At very low fields, H < H an i Sl and therefore at temperatures very near T c (i.e. \T — T c \ < 
Hanis /\dH^ 2 /dT I), the hexagonal oscillations of disappear. The AFM order parameter becomes locked 
to the lattice by the in-plane anisotropy energy as H — > 0. For H <g; H an i S , the Zeeman energy is a pertur- 
bation and leads to small deviations of rh s from the minima of the anisotropy energy. Assuming that rh s is 
oriented at = 7r/6 in zero field (corresponding to a minimum along the a* axis), the perturbed orientation 
of rh s at finite H is, 



G 



■d = ir/6- 



1 



H 



sin(2i? ff - 7r/3) . 



36 y H anis t 

This deviation leads to an orthorhombic modulation of the AFM order parameter, 



Sm s 
m 



3 /£/„ 



1 

1 - - 



H 



3 Vff a 



cos 2 (i?ff -tt/6) 



H « Bar. 



(23) 



(24) 



For# c2 (T)<if 
temperature, T c+ , 

ff c2 (T) = 5+ 



and — 7r/6 7^ ±7r/2 the upper critical field line bends towards the zero-field transition 



(r+ - t) - - 



1 + /V& 



ft 



U, 



exch 



S + {T+-T) 



COS 2 ($ H - 7I-/6) 



(25) 



In conclusion, I have developed a Ginzburg-Landau theory for the upper critical field anisotropy of UPt3 
based on a superconducting order parameter belonging to a two-dimensional representation coupled to an 
in-plane AFM order parameter. The key terms responsible for the hexagonal anisotropy of H c2 near T c are 
the in-plane anisotropy and Zeeman energies of the AFM phase. These terms lead to hexagonal modulations 
of the AFM order parameter, which generate hexagonal anisotropy in H C 2(T) through the SBF coupling to 
superconductivity. The GL theory accounts for the basic features of the observed oscillations, including the 
sign reversal at T w T* . 
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Fig. 1 Anisotropy of H^(T) vs. T/T c . The solid (dashed) curve is H^(T) for ■& H = [H\\a] (# H = tt/6 
[J?||S*]). The cross-over occurs in the vicinity of T* ~ 0.8T C . The magnitude of the anisotropy and the 
width of the transition region are amplified by the large value of Uanis/Uexch = 0.2 (chosen to show the 
anisotropy clearly in the graph). The other parameters are S~ / S + = 1.5, ertio = 0.03, (e|| +e±)mo = 0.2, 
H a nis/S + T c = 0.04. Below H anis (dotted line) the curves merge at T c (sketch). 
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Fig. 2 Anisotropy of H~^(T) vs. $h- The angular dependence of H C 2 shows the 180° relative phase of 
the oscillations for T > T* and T < T*. The parameters are the same as those in Fig. 1. The scalloped 
shape of the oscillations arises from the deviation of i?h — $ from 7r/2. 



9 



